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In this paper, a generalization of standard spin fluctuation theory is considered 
which takes into account orbital degeneracy effects which are critical for describing 
f electrons. This theory leads to an instability for a superconducting pair state 
which obeys Hunds rules. Such a state has L=5, S=l, and J=4. The degeneracy 
of this state is broken by crystalline effects, and realistic calculations for UPt% find 
a resultant pair state with symmetry, consistent with current experimental 
constraints. 
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From the beginning of theoretical work on heavy fermion superconductors, it has been 
realized that there are strong connections between these metals and superfluid 3 He. 1 This 
has led many theorists, including the author, to apply standard spin fluctuation theories 
which were developed for 3 i7e to the heavy fermion problem. So far, the results have 
been mixed. On the plus side, such theories give non s-wave pairing states, and the 
evidence in most cases is that the heavy fermion superconductors are non s-wave. On 
the minus side, the actual group representation these theories predict for UPt^, the best 
studied of the heavy fermion superconductors, has so far not matched what we think the 
experimental data are telling us. Available data point to the pair state having (£?2u) 
symmetry. 2 This state is an odd parity two dimensional group representation with line 
and point nodes, and invariably is suppressed in the spin fluctuation calculations. 3 There 
are further qualitative problems with these theories. In Table 1, a list of the seven known 
heavy fermion superconductors is shown. There are two obvious facts about this table. 
Six of the seven superconductors are uranium alloys. Second, all of the superconductors 
either have two formula units per cell, a point first remarked on by Anderson, 4 or have 
a magnetic/structural phase transition at a temperature above T c so that there are two 
formula units per cell. Another interesting point is that the magnetic susceptibilities of the 
two heavy fermion superconductors UPt$ and UPd^Al?, look almost identical to that of 
PrNi§, a localized / 2 system. Moreover, the magnetic susceptibility observed in URU2S12 
can be easily explained by an f 2 configuration. The above facts suggest that some on-site 
interaction is playing a fundamental role in heavy fermion superconductivity, since such 
an interaction could (1) differentiate between Ce and U ions and (2) depend on having 
two formula units per cell due to having in phase or out of phase relations between the 
order parameters on the two sites. 4 Standard antiferromagnetic spin fluctuation models, 
based as they are on having an attractive interaction between near neighbor sites, do not 
directly address these points. 
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The above issues indicate a need to go back and look more closely at the actual con- 
nection between the heavy fermion and 3 ife cases. The effective potential of two bare f 
electrons on a uranium site looks very similar to the direct interaction potential for ^He. In 
particular, the potential is strongly repulsive at short distances due to the direct Coulomb 
interaction, is attractive at intermediate distances (of order 3 a.u.) due to the uranium ion 
core, and decays to zero at large distances due to the exponential decay of the f electron 
wavefunction. The ground state of this potential is well known to have a symmetry of 3 #4 
(S=l, L=5, J=4) as this state minimizes the Coulomb repulsion. This represents a qual- 
itative difference between this case and that of 3 He. In the latter case, one has maximal 
S also, but as there is no orbital dependence to the bare interaction, the L state is fixed 
by the Landau parameters which are difficult to calculate. In the current case, though, 
the orbital dependence of the interaction automatically fixes the L state, with the J state 
being fixed by the strong spin-orbit interaction. 

To understand this problem further, it is useful to review the multiplet structure for an 
f 2 ion. The energy terms for the multiplets are best expressed using Racah parameters (LS 
scheme). Every configuration has an energy Eq which represents the Coulomb repulsion 
U (Eq is equal to the L = Coulomb multipole integral F° plus a combination of F 2 , 
F 4 , and F 6 terms). The splitting between singlet and triplet spin states is controlled by 
the parameter E\ (a combination of F 2 , F 4 , and F^ terms), with the 3 triplets having a 
coefficient of and 3 of the 4 singlets have a coefficient of 2 (the singlet 1 5'o is t ne highest 
energy state with a coefficient of 9). Note the similarity to the paramagnon model for 
3 -He where the coefficients are the same (0 and 2), but with the important difference that 
the splitting in the current case is not determined by the F° (charge fluctuation) term as 
in the single orbital Hubbard model used for 3 ife but by the L > (shape fluctuation) 
terms. Moreover, the degeneracy of the 3 triplets is lifted by an orbital splitting term, 
F3, which is another combination of L > terms. The lowest energy state is 3 if with an 
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energy Eq — 9E% with the next highest state being 3 F with an energy Eq. The 3 F energy 
thus represents the energy of two uncorrelated f electrons and thus sets the zero of energy 
(analogous to the 3 .He problem, where the energy of the triplet sets the zero of energy). 5 
Unlike the 3 iiZe paramagnon problem, the state has an energy lower than the energy 
zero, i.e., the interaction is already attractive at the bare interaction level. This is only 
true for a uranium (/ 2 ) ion; for a cerium ion, the zero of energy is set by the energy of 
the f l configuration, and thus the bare interaction is repulsive for all f 2 states. The bare 
binding energy is equal to the excitation energy from the 3 i?4 ground state to the 3 F2 
state. This transition has actually been seen in UPtj, by high energy neutron scattering 6 
and has an energy of 0.373 eV. 

It is a non-trivial problem to take this bare binding energy and convert it into an 
effective binding energy to use between two f electron quasiparticle states at the Fermi 
surface. To do so, collective effects must be included. First, one takes the bare interaction 
vertex to be the antisymmetrized Coulomb interaction (direct minus exchange term). For 
the s electron case, one can easily show that doing a diagram sum using this bare vertex 
leads to the standard paramagnon results. 7 The f electron case is more complicated due 
to the presence of 4 interaction parameters (Eq, E\, E2, -E3) and 14 orbitals. If one 
only keeps the £0 term, the diagram series can be analytically summed. The result is 
Eq/(1 — Eqxo)/(1 + 13-EoXo) + EqXo/(1 ~ E OXo)- This has some important implications, 
in that the effective repulsion is reduced compared to the bare Eq as long as xo is not 
too close to being equal to 1/Eq (note the difference again to the s electron case, where 
the bare repulsion is always reduced in the density channel 7 ). One can speculate that 
the divergence for xo = represents a localization instability. If one keeps only the 

E\ term, the diagram series can also be summed. For the triplet states, one obtains 
-ll^fxoA 1 - 81 ^iXo)/( 1 ~ 4E iXq) + 2-EiXo/(! - ^ e \Xq) which has some similarities 
to the s electron case. In particular, there is an induced attraction for the triplet states. 
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Note the divergence for 9E\xo = 1 which plays the role of a magnetic instability. The 
analogous series has also been done for the l I state, where the induced interaction is 
(4E 1 + 13E 2 X ~ 1264x0 " ^Ef X 3 )/(l - 81E 2 x 2 )/(l - 4E 2 X 2 ) ~ 2E V This also has 
similarities to the s electron case, and one finds for singlet states an increased repulsion. 
An analytic expression for the general series has not been obtained due to the complicated 
coefficients associated with the E<l and E% terms. Instead, one can reduce the vertex 
equations to a matrix equation which can be solved on a computer. This has been done 
for the 3 if, 3 F, 3 P, and ^1 states. Complete orbital degeneracy has been assumed and no 
spin-orbit effects have been included. 8 

The results are summarized in Figure 1, where the various effective interactions are 
plotted versus xo along with the zero of energy for the f 1 case (0) and f 2 case (Eq). The 
values of Ei were gotten from Goldschmidt^ (these values give an term of 1.83 eV, 
consistent with spectroscopic data in heavy fermion uranium compounds). As one can 
see, the triplets become increasingly attractive and the singlet increasingly repulsive as xo 
increases with a divergence for (Eq + 9Ei)xo = 1. Xo is difficult to estimate since spin-orbit 
and anisotropy effects play a major role. 10 For illustrative purposes, we assume a "Stoner" 
renormalization of 4 as seen in 3 He. For this value of xo (0.137), the 3 H energy is -2.3 eV 
relative to the f 2 zero of energy. Even for the f 1 case, there is still an (induced) instability 
for 3 if, so that pairing is indeed possible for cerium ions, although less likely. 11 

One can now estimate the effective pairing matrix element by realizing that the quasi- 
particle renormalization in the heavy fermion case is mostly frequency dependent in nature. 1 
This would then act to renormalize the induced interaction discussed above by a factor 
of Z 2 since each of the four external lines in the vertex is renormalized by Z 1 / 2 (only Z 
of the bare f electron is in the quasiparticle pole). Z _1 is equal to the mass renormaliza- 
tion factor, known from deHaas-vanAlphen measurments to be about 16 in UPt%} 3 This 
renormalizes the 3 H matrix element of 2.3 eV to about 100 K. 14 This value will be further 
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reduced when projecting onto pair states at the Fermi energy which have the symmetry 
of a particular group representation. Below, this projection factor is shown to be about 
1/8, so the final value is 12.5 K. Since the renormalized (quasiparticle) Fermi energy, Ep, 
is about 60 K in UPt% (specific heat 7, neutron scattering linewidth), the pairing coupling 
constant, NV^Z^, is about 0.21 (where N is the renormalized density of states and V4 is 
the interaction potential in the channel). With a cut-off of order Ep, this gives a 
BCS estimate for T c of 0.6 K. The agreement with experiment is somewhat fortuitous, of 
course, but the point is that the effective coupling constant is at least of the right order of 
magnitude. 

The actual symmetry of the gap is found by constructing the quasiparticle pair state 
\k,—k > using relativistic band structure wavefunctions and projecting this onto J=4 (this 
is a specific example of a general procedure advocated in the past ). The degeneracy 
of the J=4 state is broken due to lattice effects which should be well described by the 
momentum dependence of the band structure wavefunctions (although these wavefunctions 
fail to describe the frequency dependence of the quasiparticle states, they give a Fermi 
surface shape in good agreement with experimental data, indicating that their momentum 
dependence is reliable). For hexagonal UPt^, the 18 fold degeneracy of J =4 in the isotropic 
case (2J+1 times the number of f sites in the unit cell, which is two) will be broken into 
3 singlets (i~i, i~3, and i^) and 3 doublets (two J5, one Ig), with each representation 
occuring twice (+ (even parity) representations have the two atoms in phase, and - (odd) 
have the two atoms out of phase). In Table 2, these states are given in terms of pairs 
of single particle J=5/2 f states. The group transformation properties of these states are 
listed in Appel and Hertel. 15 

For each k point, there are four degenerate states available to construct \k,—k > from. 16 
The singlet (even parity) combination is (\k,Tk > —\PTk,Pk >)/2 (denoted do) and the 
triplet (odd parity) combinations are \k,Pk > (—d x + id y ), \PTk,Tk > (d x + id y ), and 
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(\k,Tk > + \PTk, Pk >)/2 (d z ), where P is the parity operator, and T the time reversal 
one. The odd parity combinations define a "d" vector which lives in a pseudo-spin space. 

The resulting pairing matrix element for this model is then 

< k', -k'\H eff \k, -k > P = (V A Z 2 )Al^ j, A^ j (1) 

where P represents the projection with A^ being the coefficient of the expansion of 
\k, —k > which has J=4 with the symmetry of the group representation r and basis v (for 
a two-dimensional representation), and j is the index of the d vector (0 for even, x,y,z for 
odd). Since the matrix element is separable in k and k', it is trivial to write down the 
appropriate BCS coupling constant 

X = (NV,Z 2 )J2<\A r k UJ \ 2 > k (2) 

j 

where N is the density of states, <>& is an average over a narrow energy shell about the 
Fermi energy, and j runs over for the even parity case and x,y,z for the odd parity case. 

The J=5/2 part of the band structure wavefunctions can be written as \k >= a ^lA* ~>i 
where \i runs from -5/2 to 5/2, % is the site index (1,2), and n is the band index (band 
calculations predict that five f bands contribute to the Fermi surface of UPt^; such a sur- 
face is in good agreement with deHaas-vanAlphen data 13 ). Thus, the A coefficients can 
be written as ^f°^ fc with k denoting either k or PTk and -k denoting Pk or Tk, with 
the appropriate linear combinations being those which match the basis states in Table 2 
and have the correct parity form (d,Q for even and d x ,dy,d z for odd). The average in Eq. 
(2) was done by constructing a regular grid of 561 k points in the irreducible wedge (1/24) 
of the Brillouin zone and keeping those nk states which are within 1 mRy of the Fermi 
energy (182 nk points for the current case). One note is that the coefficients of the I5 basis 
function in Table 2 are arbitary (subject to normalization). These are determined by a 
variational principle, i.e. one finds the set of coefficients that gives the maximum coupling 
constant. 
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In Table 3, the results of this calculation are given. The odd parity states have larger 
coupling constants since there are three terms contributing instead of the one term for the 
even parity case. This is of interest since the odd parity states only exist because of the 
presence of two f atoms in the primitive cell, which, as mentioned in the introduction, all 
heavy fermion superconductors have. The largest coupling constant occurs for a state of 

symmetry. This state is an odd parity two dimensional group representation. It has 
point nodes along the c axis and a line of nodes in the k z = ir/c zone face. It is interesting 
to note that although only the d z component of the gap function vanishes on the k z = 
zone face as expected based on group theory arguments, 17 all three d vector components 
vanish on the k z = ir/c zone face, proving a counterexample to the argument in those 
papers that a line node gap function is not possible for odd parity states. Although the 
actual form of the gap function in the current case is extremely complicated since the 
a™!? are strong functions of momentum, this state (1) is from a two dimensional group 
representation and can thus explain the unusual phase diagram seen for UPt^, (2) has the 
correct nodal structure to explain various thermodynamic data of UPt%, and (3) is an odd 
parity state with the largest possible moment projection onto the basal plane for a two- 
dimensional group representation (Mj = ±1) which is necessary to explain the observed 
directional anisotropy of the upper critical field. 18 It should be remarked, though, that 
the states and r± have coupling constants close to that of JjT and the ordering of 
the coupling constants will thus be sensitive to the cut-off of the energy shell used in the 
averaging in Eq. (2). The values tabulated in Table 3 should be multiplied by the quantity 
NV^Z^ to convert to an actual coupling constant, and as discussed above, the resulting 
coupling constant for is of the right order to explain the observed value of T c . Similar 
calculations have also been done for J=2 ( 3 F) and J=0 ( 3 -P). For J=2, the largest coupling 
constant also has symmetry (its value modulo V<i is 0.85 of the J=4 one). For J=0, 
the largest coupling constant has (s-wave) symmetry. Its value modulo Vq is a factor 
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of four larger than J = 4, so it is reassuring to find a repulsive Vo over a wide range of 
Figure 1 (in the JJ coupling scheme, Vi and Vq are repulsive). 

In conclusion, an orbital degenerate generalization of the 3 -He paramagnon model has 
been applied to f electrons and yields a superconducting pair state which satisfies Hunds 
rules (L=5, S=l, J=4). The degeneracy of this state is lifted by crystalline effects. Realistic 
calculations for the case of UPt% give a pair state with r§ symmetry which is consistent 
with experimental data with a reasonable estimate for T c . The theory also explains the 
preference for heavy fermion superconductors to be uranium alloys, and also the role that 
the crystal structure (two formula units per unit cell) plays in the pairing. 

This work was supported by U.S. Department of Energy, Office of Basic Energy Sci- 
ences, under Contract No. W-31-109-ENG-38. The author would like to acknowledge 
the hospitality of the Physics Dept., Uppsala University, where some of this work was 
completed, and to thank Kathyrn Levin for many helpful discussions. 
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Table 1. List of known heavy fermion superconductors with the number of formula units 
per unit cell. In parenthesis is the nature of the low temperature distorted phase in 
the single formula unit case (QP - quadrupolar, AF - antiferromagnetic, ? - not fully 
determined) and the resulting number of formula units. 



Case 


Formula Units 


UPt 3 


2 


UBe 13 


2 


U 2 PtC 2 


2 


URu 2 Si 2 


1 (QP/AF - 2) 


UPd 2 Al 3 


1 (AF - 2) 


UNi 2 Al 3 


1 (AF - 2) 


CeCu 2 Si 2 


1 (? - 2 ?) 
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Table 2. Hexagonal basis functions for J=4. The forms listed in this table should be 
(a) antisymmetrized > \v > —\v > |/z >) and (b) symmetrized (+ representation) 
or antisymmetrized (- representation) with respect to site before use. For A, a and f3 
are variational coefficients such that the sum of their squares is equal to one, and this 
representation occurs twice (a, (3 and (3,— a). Note that both A and A are doublets 
obtained by replacing > by | — \i >. 



Rep Basis Function 

r 5 a\5/2 > |3/2 > +/3(0.8018|5/2 > | - 1/2 > +0.5976|3/2 > 1 1/2 >) 

r 3 0.7071|5/2 > 1 1/2 > +0.7071| - 5/2 > | - 1/2 > 

A 0.7071|5/2 > |l/2 > -0.7071| - 5/2 > | - 1/2 > 

r 6 0.5345|5/2 > | - 3/2 > +0.8452|3/2 > | - 1/2 > 

A 0.2673(|5/2 > | - 5/2 > +3|3/2 > | - 3/2 > +2|l/2 > | - 1/2 >) 
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Table 3. Coupling constants for J=4 for UPt^. These are normalized to the coupling 
constant for the J=0, (s-wave) state and should be multiplied by this quantity (0.495, 
which is the square of the ratio of the J=5/2 f to the total density of states) and the 
quantity NV4Z 2 to convert to real coupling constants. Even parity corresponds to + 
representation and odd parity to -. 



Rep 


even 


odd 


A 


0.139 


0.148 


A 


0.059 


0.203 


A 


0.048 


0.129 


A 


0.027 


0.242 


A 


0.036 


0.253 


A 


0.153 


0.229 
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FIGURE CAPTIONS 



1. Effective interaction (eV) for A H, 6 F, 6 P, and I versus xo f° r parameters appro- 
priate to a U ion 9 (E = 1225 meV, E\ = 470.3 meV, E 2 = 1.923 meV, E 3 = 43.28 
meV). The zeros of energy for the f 1 and / 2 cases are marked by the dashed lines. 



